Spin polarized current in a junction of zigzag carbon nanotube 
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We investigated spin-resolved electronic transport through a junction composed of a nonmagnetic 
metal electrode and a zigzag carbon nanotube (ZCNT) by means of self-consistent Green's function 
method in the tight binding approximation and the unrestricted Hartree-Fock approximation. Our 
results show that the electric current can be spin-polarized if the coupling of the junction is weak. 
Further calculations on spin-spin correlation and local density of states reveal the existence of 
magnetic edge states in ZCNTs, which is responsible for the observed spin-polarized current and 
can be controlled by applying a gate voltage. We also studied the influence of the nearest-neighbor 
Coulomb interaction and the junction coupling strength on the spin-polarization of the current. 
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I. INTRODUCTION 

Carbon nanotubes (CNTs) have been the subject of 
a great amount of theoretical and experimental works 
during the past two decades, because of their rich and 
fascinating properties and potential applications in nan- 
otechnology. Electronic structures and transport proper- 
ties of infinite-length CNTs have been extensively stud- 
ied and well understood. A finite-length CNT, on the 
other hand, is usually treated only as a simplified one- 
dimensional quantum wire in transport theories)^ or tun- 
neling experiments (^i^i^i^ii^i^. with its crystal fine struc- 
tures largely ignored. In fact, the hexagonal lattice struc- 
ture of a CNT, inherited from the graphene lattice, may 
give rise to peculiar properties near its edges. One ex- 
ample is a recent work by Son et at, where a zigzag 
graphene nanoribbon is predicted to be half-metallic 
when an strong-enough external electric field is applied 
transversely, due to its spin-polarized edge states Sim- 
ilarly the ground state of a ZCNT also shows antiferro- 
magnetic order as a result of the Coulomb interaction 
and correlation of electrons in the hexagonal crystal lat- 
tice. This is explained by two mathematical theorems 
for the Hubbard model defined on a bipartite lattice 
with equal numbers of sublattice sites, i.e., the ground 
state is a spin singlet^^ and the antiferromagnetic cor- 
relation is always dominant over other correlations^^ 
Some literatures have investigated the band structures 
and antiferromagnetic orders (or spin density waves) in 
graphene ribbons or CNTs by the extended Hubbard 

Model3iiiil5a6il7il8il9^i21 ^Qj.]^g j^j^yg ^ggj^ jQj^g^ 

however, on the edge states of CNTs and their effects on 
transport. 

In fact, not only zigzag graphene nanoribbons but 
also finite-length ZCNTs have spin-polarized edge states 
due to the Coulomb interaction and their bipartite lat- 
tice structuresj^ The spin polarized edge states existing 
in finite or semi-infinite ZCNTs can play an important 
role in their spin transport properties. Particularly, in a 
junction composed of a semi-infinite nonmagnetic metal 



(NM) electrode and a semi-infinite ZCNT, an edge state 
may exist in the ZCNT near the interface. The edge state 
is spin polarized, with oscillations upon different sublat- 
tices, and forms a spin-dependent scattering potential. 
It will act as a spin filter when electrons tunnel through 
the interface to induce a spin-polarized current. In this 
paper, we use the extended Hubbard model to investi- 
gate the spin-dependent transport of the junction of the 
ZCNT to the NM lead. 

The paper is arranged as follows. In Sec. II we intro- 
duce the model Hamiltonian and describe the general 
formalism of our self-consistent method. In Sec. HI, the 
main results and discussions are presented. Sec. IV is a 
brief summary, and finally, an appendix is attached to 
discuss the validity of the self-consistent method we used 
in this paper. 



II. GENERAL FORMALISM 



Model Hamiltonian 



We consider a conjunction of a semi-infinite tubal lead 
and a semi- infinite ZCNT, each of which extends to an 
electron reservoir at infinity, as shown in Fig. [l^a). The 
ZCNT part is the standard hexagonal-crystal tube, and 
the NM part is assumed to be a square-lattice tube, which 
is connected to the terminal layer (layer- 1 in Fig. [IJa)) 
of the ZCNT. In the tight-binding approximation, the 
model Hamiltonian for the whole system reads 
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Here, H^j^j describes the semi-infinite NM lead, and 
df^{di„) is the creation (annihilation) operator of an elec- 
tron at the site in the NM part with spin-a (cr =|, |). 

denotes a pair of nearest-neighboring sites i and j 
in the lattice, and represents the hopping integral be- 
tween them, and Vl is the electrostatic potential in the 
NM part. Hl^^ represents the tunneling process between 
the NM lead and the ZCNT. In Fig. [^a), tf^t ^ con- 
nects the interface layers of the NM lead and the ZCNT. 
Hnt describes the 7r-orbital electrons of the zigzag nan- 
otube. We adopt the extended Hubbard model including 
the on-site Coulomb interaction U and nearest-neighbor 
interaction y ^is, 17^23 g+ (q^) jg the creation (annihila- 
tion) operator of an electron at the i*'* site in the ZCNT 
with spin-CT, = cf^Ci^, and t is the hopping integral, 
which is about 2.7eV for carbon nanotubes,'^'* taking the 
unit of energy in this paper. Vnt is the electrostatic 
potential in the ZCNT. 

To solve the interacting problem, we perform the un- 
restricted Hartree-Fock approximation, and reduce Hj^t 
to a mean- field one. 



h-NT = - E + ^ (c/CTC,,cr))c+ CjCT + VNT^riia 

+ U^iirit^) n^i + {n,{) n,^) + (n,) Uj, (3) 

where the expectation values of physical quantities (■ ■ ■) 
need to be calculated self-consistently. We use the zero- 
temperature Green's function method to establish a set 
of mean field equations, and calculate these expectation 
values for the system self-consistently. The technique 
of principal layers is used to solve the Green's function 
of a semi-infinite lead^^i^ This technique, limited to 
treat periodic structures, is suitable for the NM lead, 
but not for the ZCNT part, because the inhomogeneity 
of (c^o-Cicr) a-iid {ilia) near the interface of the ZCNT will 
make the mean-field Hamiltonian Hnt generally not a 
periodic one and thus we cannot apply this technique di- 
rectly. But we notice that the expectation values deep 
inside the ZCNT part are nearly periodic and very close 
to the corresponding values in an infinite ZCNT. So we 
divide the ZCNT into two part: the central part with N 
units (see Fig. [TJa) for the definition of 'unit' in this pa- 
per) including the interface, and the semi-infinite right 
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FIG. 1: (a) Schematic of a NM-ZCNT junction, the dashed 
rectangles denote the unit cells of ZCNT, and the dashed 
lines denote the layers of ZCNT. (b) Schematic diagram of 
the electrostatic potential and electrochemical potential. 



part. We assume the expectation values in the right 
part are periodic, in this way the Green's function of 
the right part can now be solved with the technique of 
principal layers and the whole system can be solved self- 
consistently. Thus although the whole system in Fig. [1] 
has only two compositions, in practice we treat it in three 
parts: the NM part on the left naturally, N units of the 
ZCNT including the interface at the center, and the rest 
of the ZCNT on the right. The assumption that the ex- 
pectation values {cj^Cicr) and {riia) are periodic in the 
right part, and equal to those in an otherwise integral 
ZCNT, is justified when N is so large that the expec- 
tation value of any quantity in the ZCNT changes little 
with increased N . In Appendix we will discuss this is- 
sue in more details. Based on the above analysis, the 
Hamiltonian for the ZCNT can be rewritten as 
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where hfjrp and hf^rp describe the central part and the 
semi- infinite right part, respectively, — — ^ (t + 

V (^Cj^Cicr'))cl^Cjcr rcprcscnts the connection between the 
central part hf^rp and the right part hf^rp. In this way, 
hpfrp can be treated by the technique of principal layers. 

In order to reduce the influence of band structure of 
the nonmagnetic metal lead on the transport property, we 
take tL = 4:t, which is equivalent to the so-called "wide- 
band limit" ^ And in order to preserve the edge state of 
the ZCNT, we take tf^^. — O.lt, which is a comparatively 
small value implying weak coupling between the NM lead 
and the ZCNT, or a high potential barrier between them. 
This can be realized in experiments by fabricating an 
oxide layer at the interface In equilibrium, the Fermi 
level of the whole system Ef is identical everywhere as 
shown in Fig. [^b). Following the definitions in Ref. 
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|2a |. the electrochemical potential (ECP) in the ZCNT 
is defined by /xjvt = ~ VnTj and the electrochemical 
potential in the NM part is defined by fi^ = Ef — Vl- 
In this paper, Ef, Vl, and /i^ are fixed, Vnt is assumed 
to be adjustable by means of a gate voltage. Adjusting 
fj-NT may control the electron density He in the ZCNT. 
We will fix the NM part at half-filling by tuning Vl (or 
/ii), so that the transport modes in the left lead around 
the Fermi level can match those in the ZCNT, to facilitate 
the electronic transport 



B. Method 

We come to establish the mean field equation by means 
of the self-consistent Green's function method, and then 
to use the recursive method to find the numerical solu- 
tion. The retarded Green's function is defined by 

GL,,,, {t;t') = -te{t-t'){{c,^ it) , c+ it')}), (5) 

whose Fourier transform can be calculated by the Dyson's 
Equation's 
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where rj is an infinitesimal constant, and h^j. is the 
mean-field Hamiltonian of the central part as defined in 
the last subsection. El = H^^g^H^^ is the self-energy 
due to the left lead, where H^^^ is the Hermitian conju- 
gate of and g}^ = [E + ir] — Hj^j^)~^ is the retarded 
Green's function of the left lead. The self-energy due to 
the right lead has the similar definition S/j = hf^g^h^^ 
and g^i^ = {E + irj — hpjrp)^^. In the present approach, 
the system is defined on a lattice and the Hamiltonian 
hfjrp as well as the self energy can be re-expressed in the 
form of square matrix. The dimensionality of the matrix 
is twice of the lattice number of the cental part of car- 
bon nanotube because the spin degree of freedom. For 
a finite number of lattice sites, the Green's function can 
be calculated numerically once the initial values are as- 
signed to the expected quantities in Eq.dS]) such as the 
local densities of charge and spin. Once we have this 
Green's function, the expectation value of any physical 
quantity can be obtain as follows 
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where G[* j^, {E) is the complex conjugate of G,[^ j^, {E). 
Then the local charge and spin density are given by 



and 



(8) 
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respectively. 

The recursive method is applied to obtain the stable 
solution for the problem. The stable solutions of the 
Green's functions are applied to calculate the relevant 
physical quantities. The spin polarization (SP) of the 
conductance is defined by 



Sp = 



Gt — G 



Ct + Gj 
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where G^ is the conductance in the spin-cr channel (cr =t 
, I), which is given by the Landauer's formula, 
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'Tr[rLcrGTlJcrG^^]E=Ef , 
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where Ta^ = i{^Aa — ^^ao-)' ^^'^ ^^<^ represents the self- 
energy of lead A (= L, R) in the spin-cr subspace, G''+ is 
the Hermitian conjugate of G^ . 

For concreteness, we describe the overall procedure 
of the present self-consistent calculation, which is illus- 
trated in a fiow chart in Fig. [5] and elucidated as follows: 

Step 1: Calculate the self-energy due to the left lead. 
Because there's no Coulomb interaction term in the 
Hamiltonian of the left lead, we can solve this part with- 
out recursion. An arbitrary value can be assigned to 
Vnm to complete the Hamiltonian Hf^j^j, then g£(-E) and 
Si(i?) can be calculated by the technique of principal 
laye r^^i^^ . It's found that the left lead is half-filled when 
the Fermi energy Ef = Vl, as the left lead is a simple 
square lattice. The Fermi level Ef is identical for the 
whole system because we are studying in the linear re- 
sponse regime, and it is kept unchanged in the following 
steps. 

Step 2: Construct a trial Hamiltonian h^T- An ini- 
tial charge and spin distribution is assigned in the zigzag 
nanotube. The initial change density promises that the 
number of electrons per site is about 1, which represents 
the half-filled case of the zigzag nanotube. Spin density 
wave and charge density wave are constructed on the 
zigzag nanotube at the same time, because many litera- 
tures have referred to the coexistence of the spin density 
wave and charge density wave in carbon nanotubes or 
ribbon£^^'^-i2£i2ij and this is also a general property of 
the half-filled extended Hubbard modeP^'"^^. The hop- 
ping term (c^Qcr) is assigned to be 0, because we have 
no knowledge about it so far. Remember {^c^^Cia) and 
(uicr) in the right unit of the central part is also used to 
determine the Hamiltonian and h^j,. The electro- 
static potential Vnt can be assigned the same value with 
Ef at first, because all cases are nearly half-filled. Vnt 
will be changed in the following iterations. Although Vnt 
is fixed by the gate voltage in experiments, we cannot use 
it to be the restrictive condition in the calculation, oth- 
erwise the iterations won't converge. We use the local 
charge density Uf, in the right lead to be the restrictive 
condition, which can lead to convergence. Because rig 
has a one-to-one correspondance to Vnt, Vnt will be 
determined along with the procedure of convergence. 
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timated from some organic molecules or by fitting the 
exciton energies of nanotubes. People also try to use 
the first principle calculations to obtain the values of U 
and V . But as the authors of the Ref. [20*1 pointed out, 
the values of these parameters depend on the choice of 
the exchange-correlation functional used within the cal- 
culation of density-functional theory. For example, U is 
chosen at QM ~ 2t when one studied the problem of edge 
states'^^. Here we adopt this result to discuss the effect 
of U and V in several cases. 



III. RESULTS AND DISCUSSION 




FIG. 2: The flow chart of the self-consistent calculation. 



Step 3: Calculate the Green's function G"'(i?) in 
Eq.®. The self-energy Yjb{E) can be obtained from the 
Hamiltonian hf^^ and h^rp. 

Step 4: Calculate the expectation values of the physical 
quantities (cj'^Cicr') and {uicr) from Eq.© and ([5]), and 
find a new Fermi energy Ey The restrictive condition 
for finding E'j: is the local charge density Ue in the right 
lead is a certain value. After that we have to draw the 
Fermi energy back to the initial value Ef by means of 
tuning VffT- So we substitute Vnt ~ {E'f — Ef) for V^t- 

Step 5: Judge the convergence. If the expectation val- 
ues of the physical quantities {c^^Ci^) and {rii^) has con- 
verged, we can end the iterations and get the stable solu- 
tions. Otherwise we have to use Vnt, {cj^Ci^) and {ni„) 
to update the Hamiltonian h^T, then turn to step 3 to 
begin a new iteration. 

Step 6: After convergence, we use the the stable solu- 
tions of the Green's functions to calculate the spin polar- 
ization of conductance by Eq. pUj) and pT|) . 

When we calculate the expectation values from Eq.©, 
for the sake of precision, we employ the contour integra- 
tion in the complex energy plane^, instead of direct inte- 
gration along the real energy axis, because the retarded 
Green's function can be analytically continued into the 
complex energy plane. On the integration path, the com- 
plex energy has a large imaginary part, so the infinitesi- 
mal constant r] in Eq.® becomes negligible, even can be 
assigned 0. But when we calculate Eq. pTj) . the constant 
•q becomes indispensable. It should be assigned a value 
as small as possible, so that the conductances are not 
influenced by it. 

On the magnitudes of U and V , U was adopted in 
a range of 2t ^ At and V k, t for graphene, ribbons 
or nanotubes in some literature a^'^'^^i^^ . They were es- 



A. The case of V = 

In this section, we first consider the case of = 
and focus on the effect of [/, since it is J7, rather than 
V , that is one of the essential factors to stabilize the 
spin-polarized edge states. We choose U — 1.2t in this 
subsection. And the ZCNT's configuration is (6,0), which 
is metallic when the Coulomb interaction is ignoredii 

Fig. [3Ka) shows the spin-resolved conductance as a 
function of the relative electrochemical potential of the 
ZCNT, defined by A/ijvT = Mjvt— Mjvt! where /i^g- 
the electrochemical potential of the ZCNT in the half- 
filled case. And Fig. ^h) shows the spin polarization of 
the conductance, calculated from the data in Fig. EJa). 
We divide the electrochemical potential into three parts, 
according to the spin polarization. In region A and C, 
the spin polarization of the conductance is zero, while 
the spin polarization is not equal to zero and even may 
change its sign in region B. In order to understand this 
phenomenon, we investigate the energy levels of the edge 
of the ZCNT, i.e., layer 1 in Fig. [T^a). 

In Fig. m the square (diamond) markers denote the 
energy level of the edge state with spin-t(|), which are 
obtained by locating the peaks of the spin-resolved local 
density of states (LDOS) in layer 1. The abscissa rig 
is the electron density deep inside the ZCNT, by our 
previously stated assumption it takes the value of the 
electron density at the right side of the central ZCNT, 
or equally in the right ZCNT part. The zero point of 
energy is chosen to be at ^nt = Matt- straight line 
describes the dependence of the relative electrochemical 
potential A^nt on the electron density Ug. The linear 
relation between them implies a constant density of states 
(DOS) and linear dispersion relations around the Fermi 
level of the ZCNT. We again divide the whole range into 
three parts according to the relative magnitudes of the 
energy levels of the edge states and the electrochemical 
potential of the ZCNT. We notice that this division is in 
accordance with Fig. [31 i.e., the energy levels of the edge 
states overlap in region A and C, but split in region B. 

In region A, Afij^x is below the energy levels of two 
edge states, which are both empty. The charge density rig 
in layer 1 is very low as shown in Fig. EKa), and the spin 
density S{i) is zero everywhere as shown in Fig. [D^b). 
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FIG. 3: (Color online) (a) Spin-resolved conductance and (b) 
Spin polarization of the conductance of the NM-ZCNT junc- 
tion, as functions of the relative electrochemical potential of 
the ZCNT. The division of the region A, B and C is according 
to the magnitude of the SP. 



FIG. 5: (Color online) (a) Electron's charge and (b) spin den- 
sity distribution along the ZCNT in the three configurations 
A, B and C. 
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FIG. 4: (Color online) The energy levels of the edge states 
and the relative ECP in ZCNT as functions of electron density 
in the ZCNT. The whole range is divided into three regions 
according to the relative magnitudes of the energy levels of 
the edge states and the ECP of ZCNT, in accordance with 
Fig. El 



Thus the charge current is not spin polarized through 
the junction in region A of Fig. [S] In region B, Afi^T 
is between the energy levels of two edge states, thus only 
one of them is occupied. This leads to a high spin density 
in layer 1 with its charge density nearly equal to that in 
other layers as shown in Fig. [Slja). The spin density 
oscillates in antiferromagnetic order in the next several 
units and decays rapidly as shown in Fig. [5l^b), which is 
the same as expected in Refs. jTllfl3 |. The spin-polarized 
edge state acts as a spin-dependent barrier, and makes 
the current through the junction spin-polarized in region 
B. In region C, /S.^mt is above the energy levels of both 
edge states, which are occupied. The charge density in 



layer 1 is very high as shown in Fig. [5ja), and the spin 
density is zero everywhere once again in Fig. El^b). So 
there's no spin-polarized current through the junction. 

Furthermore, Fig. [3l^b) shows an approximate anti- 
symmetry, but not an exact one. This fact reflects the 
partial particle-hole symmetry of the Hubbard model on 
a biparticle latticeiiiii^ As we can see in Fig. [TJ either 
the NM part or the ZCNT part is of biparticle lattice, 
but the tunneling terms between them H^^^ break the 
bipartite structure of the whole system, and the symme- 
try is broken. However because of a weak H^^^, the SP 
still exhibits the antisymmetry approximately. 

This subsection contains the main result of this paper, 
it describes three phases of the edge of the ZCNT. They 
can emerge consequently by increasing the electron den- 
sities or the electrochemical potentials of the ZCNT. As 
a result the transition between these phases can be real- 
ized by adjusting the gate voltage applied on the ZCNT 
in experiments. 



B. The case of 1/ / 

In this subsection, we investigate the effect of the 
nearest- neighbor Coulomb interaction V . We take U — 
2t, and compare the spin polarization in the cases of 
y = 0, 0.05i, O.lt, 0.2t, and 0.4<. Fig. [H shows the 
SP of the conductance as a function of the electron den- 
sity in the ZCNT, which indicates that V has two effects 
on the SP. 

(i) When V increases, the region of non-zero Sp moves 
leftward. This means that the energy levels of the edge 
states become lower from the middle of the whole band 
as V increases, which can be understood from the band 
structure of the zigzag nanoribbons in Ref.[l3| because 
the energy levels of the edge states of ZCNT can be ob- 
tained from the flat bands of the zigzag nanoribbons by 
means of the band- folding approach. ^i The flat bands in 
Fig. 2(c) of Ref. fioj are lowered from the middle of the 
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FIG. 6: (Color online) Spin polarization of conductance as a 
function of electron density in the ZCNT for different values 
of the nearest-neighbor Coulomb interaction strength V, with 
vertical offsets of 10% for clarity. 



whole band, and the same mechanism also happens in 
ZCNTs. Now we know it attributes to the effect of the 
nearest-neighbor Coulomb interaction V. 

(ii) When V = 0, the spin polarization of the conduc- 
tance is nearly antisymmetric. The position of the sign 
change is about at n,, — 1, i.e., the half- filling point as 
shown in Fig. [3l This property is relevant to the partial 
particle-hole symmetry of the Hamiltonian Eq. Pc)) at 
V = Ofiiii^ and the weak tunneling limit. But the anti- 
symmetry is broken as V increases, the region of S'p < 
expands leftward in Fig. [6l Thus we may conclude that in 
the cases with V — 0.2t, 0.4t, as in Fig. [51 the lack of any 
sign change in the spin polarization of the conductance 
is a result of the nearest-neighbor Coulomb interaction 
V. 



C. Effect of the coupling strength 

In this subsection, we investigate the effect of the cou- 
pling between the ZCNT and the NM part. Fig. [7] illus- 
trates the energy levels of the edge states and the relative 
ECP in the ZCNT as functions of the coupling strength 
between the ZCNT and NM electrode. The relative ECP 
is always zero because the ZCNT we study in this subsec- 
tion is always half- filled. We can see the energy levels of 
the edge states shift down when the strength of coupling 
increases. There's a sudden change around tf^^ = 0.17t, 
where the two energy levels of the edge states are both 
occupied. This is a transition from the single to dou- 
ble occupancy on the edge of the ZCNT. Because the 
edge states are spin-polarized, the transition is also ac- 
companied by a magnetic transition at the edge. The 
inset shows the spin density on the edge of the ZCNT 
as a function of the coupling strength if^^. When the 
edge states are single-occupied, the edge of the ZCNT 
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Strength of coupling between ZNT and NM (t) 

FIG. 7: (Color online) The energy levels of the edge states 
and the relative ECP in ZCNT as functions of the coupling 
strength between the ZCNT and NM electrode. The inset 
shows the spin density on the edge of the ZCNT as a function 
of the coupling strength. The parameters are U = 1.2t, V — 0. 



is magnetic. However, the magnetism of the edge disap- 
pears suddenly at tf^^ = 0.17t when the two edge states 
are both occupied, and the spin polarizations of the two 
states cancel each other. 



IV. CONCLUSION 

In summary, we have investigated the spin polarized 
transport in a junction of a nonmagnetic metal lead and 
a zigzag carbon nanotube. Due to Coulomb interaction 
of electrons, the zigzag carbon nanotube indicates strong 
antiferromagnetic correlation at half-filling, which is en- 
hanced near the edge. The effect of the magnetic edge 
states leads to spin polarization of the electronic trans- 
port. 

Compared with the half-metallic zigzag graphene 
nanoribbons in Ref. Jfl], the NM-ZCNT junction device 
we present here can be fabricated more easily because of 
the developed synthesis technique of carbon nanotubes, 
though the spin polarization is weaker than that of the 
zigzag nanoribbons, which can be 100% because of their 
half-metallic property. In both cases edge states are the 
keys in generating spin polarization, but they act in dif- 
ferent ways. The edge states of a nanoribbon carry spin- 
polarized currents, while the edge state in our nonmag- 
netic metal-zigzag carbon nanotube junction does not, 
actually it blocks the current in a spin-dependent manner 
to produce spin polarization. The spin filter function of a 
ZCNT hasn't been discovered in experiments so far. This 
is possibly because the contact between the nanotube and 
the electrode has allowed the current to enter the nan- 
otube without crossing the edge, as shown in Fig. 1 of 
Ref. ^4 . This type of contact is called "side-contact" in 
Ref. |35j . What we need is the so-called "end-contact", 
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FIG. 8: (Color online) Spin densities in the leftmost and right- 
most layer of the central part as functions of the lengh of the 
central part for U — 1.2t and 2t, V = 0. The spin density in 
the rightmost layer is multiplied by (-1). 
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Appendix: Validity of the self-consistent method 



In the Sec. n(A), in order to use the technique of the 
principal laye r^^i^^ to calculate the Green's function of 
the semi-infinite ZCNT lead, we divide the NM-ZCNT 
junction into three part: the left NM lead, the central 
ZCNT part, and the right ZCNT lead. We take the left 
N units of the ZCNT as the central part, and the rest of 
the ZCNT as the right semi-infinite lead. It is assumed 
that the right ZCNT lead has the mean-field Hamilto- 
nian hf^j. with periodic distribution of local charge and 
spin densities. This assumption becomes reasonable only 
when the length of the central part N is large enough 
such that the edge effect vanishes completely to the bulk 
properties. For a large N, the expectation of any physical 
quantity in the ZCNT lead is assumed to be independent 
of N. We investigate the spin densities in the leftmost 
and rightmost layers of the central part as functions of 
N in Fig. [8l It can be seen that the expectation values 
converge quickly as N gets larger. The spin density in the 
leftmost layer is large to form a magnetic edge state. The 
spin density in the rightmost layer, which is multiplied 
by -1 in Fig. [SJ is small, and converges to the ampli- 
tude of the spin density wave in an infinite ZCNT^^. By 
the way, we have used the self-consistent Green's func- 
tion method on the infinite nanotubes, and recovered the 
results in Ref. Ilal. 



FIG. 9: (Color online) Spin polarization of conductance as a 
function of the relative electrochemical potential of ZCNT for 
different length of the central part N=8,12, 16,20, 24. 



and the coupling between the ZCNT and the NM elec- 
trode must be weak enough to preserve the edge state of 
the ZCNT. This way the junction may act as a spin filter 
and produce spin-polarized tunneling current. 



Further more, we compare the SP of conductance for 
different length of the central part, as shown in Fig. [9] 
As N increases, most of the SP do not change except the 
region of sign reversal. The peaks of SP become lower 
as N increases. Due to the restriction of our computer's 
ability, we just take iV = 24 for the central part in this 
paper, and keep in mind the SP of conductance we calcu- 
lated may not be very accurate around the point of sign 
reversal. 
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We also make an error estimation of the self-consistent 
method. We run 10 more iterations after the condition of 
convergence has been satisfied, and calculate the relative 
standard deviation of Sp. The relative standard deviation 
are all less than 1% for dozens of random samples. 
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